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Abstract 



The chiral symmetry relation and scaling of the overlap fermions are studied 
numerically on the quenched lattices at 3 couplings with about the same phys- 
ical volume. We find that the generalized Gell-Mann-Oakes-Renner relation is 
satisfied to better than 1% down to the smallest quark mass at moa = 0.006. 
We also obtain the quark mass from the PCAC relation and the pseudoscalar 
masses. The renormalization group invariant quark mass is shown to be fairly 
independent of scale. The n and p masses at a fixed m^/mp ratio indicate 
small 0(a 2 ) corrections. It is found that the critical slowing down sets in 
abruptly at a very small quark mass close to those of the physical u and d 
quarks. 
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Recent advances in chiral fermion formulation on the lattice hold a great promise in 
implementing chiral fermions for QCD at finite lattice spacing y|. There are distinct ad- 
vantages over the previous formulations. For example, the Wilson fermion breaks chiral 
symmetry at finite lattice spacing a and, therefore, the task of extrapolating the Monte 
Carlo results to the continuum and chiral limits requires fine tuning and is often difficult. 
In this case, one could not expect low energy theorems to be reproduced at finite lattice 
spacing. Similarly, the staggered fermion can only be formulated with four flavors with an 
U(l) subgroup of the flavor non-singlet chiral symmetry. At finite lattice spacing where 
the numerical calculations are carried out, the flavor symmetry is broken. Furthermore, the 
anomalous chiral Ward identity does not hold for these fermions unless at the continuum 
limit. This makes the analysis of anomaly on the lattice rather unclear. There is no unam- 
biguous identification of the fermion zero modes with the topology of the background gauge 
field. These difficulties have rendered the studies of low energy phenomenology of QCD on 
the lattice unsettling. 

The picture has been altered dramatically with the advent of Neuberger's overlap 
fermion || which is derived from the the overlap formalism ||. All of the above-mentioned 
impediments can be avoided pending pristine numerical simulation. It is shown to have 
correct anomaly and exact chiral symmetry on the lattice and there are no order 

a artifacts ||. The overlap fermion has a compact form in four dimension and is easily 
employed to derive low energy theorems, chiral symmetry relations, and anomaly at finite 
lattice spacing. In this letter, we shall study the overlap fermion numerically. We test chiral 
symmetry via the Gell-Mann-Oakes-Renner relation at finite lattice spacing and check the 
scaling behavior of the ir and p masses. The preliminary results were reported earlier ||. 
We also obtain the quark mass from the Chiral Ward identity and verify that it is free 
from additive renormalization and the renormalization group invariant quark mass is indeed 
independent of scale. 

Neuberger's Dirac operator has the following form for the massive case M 

D(,n ) = 1 + =55 + (1 - (1) 

where e(H) = H/ \J H 2 is the matrix sign function of H which we take to be the Hermitian 
Wilson-Dirac operator, i.e. H = 75!?^. Here D w is the usual Wilson fermion operator, 
except with a negative mass parameter which corresponds to k c < k < 0.25. We take 
k = 0.19 in our calculation. The massless operator -D(O) is shown || to satisfy the Ginsparg- 
Wilson relation ||: {75, .D(O)} = D(0)^f 5 D(0). The bare mass parameter m is proportional 
to the quark mass without additive constant which we shall verify. 

There are several numerical approaches to approximate the sign function e(z) flOHTB" 



We adopt the optimal rational approximation fl2| with a ratio of polynomials of degree 12 
in the Remez algorithm. We find the error to the approximation of e(z) to be within 10 -5 
in the range [0.02, 2] of the argument z. In the range [0.0005, 0.02], the error can be as 
large as 1%. To improve the accuracy of e(H) and hence the chiral symmetry property, the 
smallest 10 to 20 eigenvalues of H 2 with eigenvalues of \H\ less than 0.04 are projected out 
for exact evaluation of the sign function from these eigenstates [0. This has the added 
benefit of reducing the number of iterations for the multi-mass conjugate gradient inversion 
of H 2 + q in the inner loop by a factor of 3.5 or so. We checked the unitarity of the matrix 
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V = j5e(H). For Vx = b, we find \x'x — tfb\ ~ 10 9 . Since V is unitary we exploit the 
identity (1 + V^)(l + V) = 2 + V + V f in order to use the conjugate gradient algorithm 
on the hermitian matrix V + instead of V^V which has a higher condition number. 
Furthermore, since [V + ^,75] = 0, one can use a chiral source, i.e. 756 = ±6 to save one 
matrix multiplication jTJJ] per iteration. In the present study, we consider three lattices: 
6 3 x 12 at p = 5.7, 8 3 x 16 at (3 = 5.85, and 10 3 x 20 at (3 = 6.0, which have about the 
same physical volume. With residuals at 10~ 7 , the inner loop takes typically ~ 200 — 250 
iterations and is almost independent of the lattice volume. The outer loop takes ~ 40 — 110 
iterations from the small to large lattice volumes. Even for topological sectors with charge 
Q 7^ 0, we find the critical slowing down is much milder than that of the Wilson fermion and 
we found no exceptional configurations. The critical slowing down sets in quite abruptly 
when m a < 0.006 for the sector with topology. This is already very close to the physical u 
and d masses. 

The lattice chiral symmetry is reflected in the generalized Gell-Mann-Oakes-Renner 
(GOR) relation 



m a / d 4 x(7i(x)n(0)) = 2(^)- 



(2) 

Since it is satisfied for each quark mass, lattice volume and spacing, configuration by config- 
uration, and for each source p,|14|,|15j, it serves as an economic test of chiral symmetry. Here 
7r(x) = V ; 75(' r /2)V ; is the pion interpolation field and the quark propagator is the external 
one [|7|,|T6| with D^ 1 = (1 — moa/2)~ 1 [D~ 1 (mo) — 1/2]. Alternatively, one can use the bilinears 
(1 — moa/2) _1 ^r(l — D/2)ip for the operators and D~ l as the propagator. We utilize this 
relation in Eq. (^) as a check of the numerical implementation of the Neuberger operator. 
We find that for the lattices we consider the GOR relation is satisfied very well (to within 
10~ 3 ) all the way down to the smallest mass m^a = 2 x 10 -4 for the Q = sector. For the 
Q 7^ sector, the presence of zero modes demands higher precision for the approximation 
of e(H). For example, we show in Fig. |1] the ratio of the right to left side of Eq. (0) for 
1.2 
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FIG. 1. Ratio of the right to left side of Eq. (g) for a configuration with topology. The symbols 
»/o indicate the case with projection of 10/20 smallest eigenmodes. 



a typical configuration with topology on the 6 3 x 12 lattice at (3 — 5.7 as a function of 
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the quark mass mod. When only 10 smallest eigenmodes of H 2 are projected, we see that 
the ratio deviates from one for small quark masses. For the smallest mass m Q a = 0.006, 
it may deviate as much as 16%. The situation is considerably improved when 20 smallest 
eigenmodes are included where the deviation is reduced to 1%. Our result with the overlap 
is appreciably better than the domain-wall fermion case when the size of the fifth dimension 
is limited to L s = 10 to 48 UTTf . In the latter case, the ratio deviates from unity by ~ 55% 
for L s = 10 and ~ 15% for L s = 48 at the quark mass m/a = 0.02 which is about 3 times 
larger than the mass of m a = 0.006 in our study where the deviation is at the 1% level. 

The average of u and d quark masses in the MS scheme at the renormalization scale \x 
is obtained from the axial Ward identity Z^d^A^ = 2Zg 1 moZpP via the ratio of the matrix 
elements 



m 



MS 



(//) = Z 5 1 (/i)m = lim 



Za{^) E^(0|V^4(x)|7r(0)) 
ZpQi) 2£a<0|P(a;)|7r(0)) 



(3) 



where = ^7^75 (t /2)ip and P = ipi^^r / 2)ip . The resumed cactus diagram of one-loop 
calculation of Za and Zp |L8| are used in Eq. (Q). It turns out that the ratio Za/Zp is very 



close to unity for /i = I /a. The renormalized quark mass a defined in Eq. ([3]) is plotted 
in Fig. against the bare mass m^a for the three lattices. 
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FIG. 2. Renormalized quark mass vs. the bare quark mass on the three lattices. 



We first observe that the renormalized quark mass does not have an additive part due to 
the lattice chiral symmetry. The linear fit including the smallest 5, 7, and 8 quark masses 
for (3 = 5.7, 5.85, and 6.0 respectively shows that the intercepts are of the order of 10~ 3 and 
are consistent with zero. This is a distinct advantage over the Wilson fermion whose quark 
mass is subject to additive renormalization which depends on the gauge configuration. From 
the slope we can determine the non- perturbative Z s l which are 0.95(5), 0.89(5) and 0.80(1) 
for (3 = 5.7, 5.85, and 6.0. These are within 14% from the perturbative calculation which are 
0.818, 0.824, and 0.829 respectively [IE]. The renormalization group invariant quark mass 



is defined as the integration constant of the evolution equation such that 



mf GI = AZ s (/x)mf 5 (/i) = AZ 5 (/i)^ 1 (/i)m . 



(4) 



Using the four- loop calculation of AZs(/J>) in the MS scheme | 2D| , we obtain the product of 
AZs(fi)Zg 1 (fj,) which should be scale invariant. The results from the 3 lattices give 1.12(1), 
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1.17(6), and 1.14(5) for (3 = 5.7,5.85, and 6.0. They are indeed independent of scale within 
errors. From the pion mass fit to be discussed later, we determine m = 6.8(5) MeV from 
the physical pion mass on the f3 = 6.0 lattice. This gives m^ GI = 7.6(6) MeV for the average 

of u and d quark masses, which in turn gives m^ 1 (p = 2GeV) = 5.5(5) MeV. The lattice 
scale is set by the r as derived from the static quark potential |2"T| . We should point out 
that this quark mass should not be taken literally, since the volume is quite small. The 
primary purpose of the present study is to verify that the quark mass extracted from the 
overlap fermion action does not have additive renormalization and that m^" is indeed scale 
invariant. 

Next we look at hadron masses. The results for the 7T, p and nucleon masses on the 
10 3 x 20 lattice at (3 = 6.0 are given in Fig. |^. They are obtained by a single exponential fit 
with covariance matrix. 
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FIG. 3. Masses of it, p and N on the 10 3 x 20 lattice at j3 = 6.0 vs. mo- The scale tq is used 
for conversion to physical units. 



We see clearly that the nucleon mass suffers from the finite volume effect when mo is 
smaller than 0.4 GeV. Although the p and tt masses appear not affected as much, there is the 
worry that the finite volume effect is present when one considers the pion mass behavior as a 
function of the quark mass mo- In the finite size scaling region where the lattice size is much 
smaller than the pion Compton wavelength (L « l/m w ), the eigenvalue distribution and 
the chiral condensate can be described by the universal function derived in the chiral random 
matrix theory |[22|| . While the chiral condensate of the overlap fermion in this region has 
been studied [p3| |, we shall concentrate on the region L > l/m n where the chiral perturbation 



analysis is expected to apply. Plotted in Fig. are the pseudoscalar meson mass squared 
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{m P a ) as a function of m^a for the three lattices for those pseudoscalar mesons whose 
Compton wavelengths are less than the respective lattice size. 
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FIG. 4. Pion mass squared as a function of the bare quark mass on the three lattices. 



We fit them with a form suggested by the quenched chiral perturbation theory [24 



2 „2 



2Am a 2 {l - Sln(2Am / A 2 )} 



ABm 2 . 



(5) 



We find that for A x in between 0.6 GeV and 1.4 GeV, the chiral log 8 for the (3 = 6.0 case 
is in the range of 0.20(2) to 042(14) which is slightly larger than what is obtained from 
the Wilson fermion |]25|| . The fit is quite stable with small x 2 {x 2 /DF < 0.1) and is fairly 
insensitive to the range of the fitted quark masses shown in Fig. || On the other hand, the 
fit for j3 = 5.85 and 5.7 lattices shows that 5 is consistent with zero (the typical value ranges 
from -0.04(21) to -0.34(48)) and the result is more sensitive to the range of quark masses 
that are fitted. We are not able to draw any definite conclusion from these inconsistent 
results, except to point out that perhaps this is due to the finite volume effect. After all, 
the physical volumes of the lattices we use are quite small (L ~ lfm). Another possibility 
is that the smaller quark mass cases may already fall into the finite size scaling region where 
the pion mass may have a different behavior than prescribed in Eq. This will be studied 
elsewhere. 

Finally, we examine the scaling of p and tt masses. Since m^a 2 and m p a are fairly linear 
in m a, we fit them with m p a = A + Bm a and m 2 a 2 = Cm a + Dm^a 2 for simplicity. 
The fits are decent with \ 2 /DF < 1 for the three lattices and the full range of quark 
masses in Fig. [|. From the fits of m p a and m 2 a 2 , we determine m for which the ratio of 
m^/rrip = 0.4, 0.5, 0.6 and plot in Fig. ^| the corresponding m p and m n in units of ^Ja as a 
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FIG. 5. Pion and rho masses in units of \J~a on the three lattices are plotted against aa 2 . They 
are given at three different m n /m p ratios. The dashed line is the fit with a constant. 



function of aa 2 , where a is the string tension. The errors on the vector and pseudoscalar 
masses are determined from interpolating the data from the neighboring quark masses. It 
is known that the overlap operator does not have 0(a) artifacts ||. It appears from Fig. |5| 
that even the 0(a 2 ) errors are small. 

To conclude, we find that when the matrix sign function t(H) in the overlap fermion is 
well approximated, the promised chiral symmetry at finite lattice spacing and the scaling of 
the renormalization group invariant quark mass and hadron masses are manifested in the 
present numerical calculation. One drawback of the overlap fermion is the large numerical 
overhead in the present algorithm. But, the stake of being able to implement chiral symmetry 
at finite lattice spacing is high. Furthermore, the unexpected feature of being able to push 
the critical slowing down to close to the physical u and d quark masses has the advantage 
of being able to study the correct chiral behavior. The nice scaling result is encouraging 
for controlling the continuum extrapolation and may afford the possibility of working at 
relatively large lattice spacings. For the moment, the study is limited to small volumes. As 
long as one can extend it to large volumes, it appears that one will be at a stage of putting 
all the systematic errors of at least the quenched approximation under control. 
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